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Abstract 

We rederive here in a simple and transparent way the master for- 
mula for the dominant part of large relativistic corrections to the 
positronium decay rate. 
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1. The strong disagreement between the experimental value of the or- 
thopositronium decay rate |IJ 

T°- p Ps = 7.0482(16) ^ (1) 

and its theoretical value which includes the order a and a 2 log(l/a) correc- 
tions [2-5] 



1 - 10.28- - -a 2 log- 
it d a. 



7.03830 fis- 1 (2) 



is a real challenge to the modern QED. For the disagreement to be resolved 
within the QED framework, the correction ~ (a/ir) 2 , which has not been 
calculated completely up to now, should enter the theoretical result @ with 
a numerical factor 250(40), which may look unreasonably large. 

Such a hope is not as unreasonable however. One class of large second- 
order corrections arises as follows ||. The large, —10.28, factor at the a /it 
correction to the decay rate (see @) means that the factor at the a/n cor- 
rection to the decay amplitude is roughly 5. Correspondingly, this correction 
squared contributes about 25(a/7r) 2 to the decay rate. Indeed, numerical 
calculations |7], |8|] give factor 28.86 at {a/irf in this contribution. 

One more class of potentially large contributions to the positronium decay 
rate is relativistic corrections. A simple argument in their favour is that the 
corresponding parameter {v/cf ~ a 2 is not suppressed, as distinct from 
that of usual second-order radiative corrections, {a/nf, by the small factor 
1/tt 2 ~ 1/10. 

This problem was addressed in Refs. [9-11]. However, the discrepancy 
between the results obtained in ||, on one hand, and in [|K], [lTJ, on an- 
other, is huge. While according to @, the relativistic correction constitutes 
(in "radiative" units, with 1/tt) 24.6{a/nf, the results of (TO, O, obtained 



in different techniques, are in a reasonable agreement between themselves: 
46(3) («/7r) 2 |10| and 41.9 (a/n) 2 [[□]]. The main source of this discrepancy 
can be traced back to the different treatment of (f/c) 2 arising in the expan- 
sion of the annihilation kernel. While the prescription of 0] is effectively (see 
their formulae (12), (18)) 



our master formula is 



{v/cf — -a 2 /4, (3) 
{v/cf — > -3a 2 /4. (4) 
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Let us mention here that this our recipe refers to the relativistic correction 
to the annihilation kernel itself; as to the phase space correction, we use in 
it the prescription (|3|). 

In view of the mentioned discrepancy, we believe that it is instructive to 
present an alternative derivation of formula (P, more obvious and transpar- 
ent than our original one. 

2. When calculating the decay amplitude, we have to integrate the anni- 
hilation kernel M(p) over the distribution of the electron and positron three- 
momenta p. We address in this note the relativistic corrections to M(p) only, 
i.e., we take as the ground-state wave function ip{p) the nonrelativistic one. 
Then the decay amplitude is 

where a = 2 /ma is the positronium Bohr radius. To lowest approximation 
in v/c the kernel M(0) is independent of those momenta and we are left with 

Thus, in the limit p — > we obtain the common prescription: the positronium 
decay rate is proportional to \ip{r = 0)| 2 . 

However, already to first order in (p/m) 2 the momentum integral 

J dpip/rnf J a ™\ )2 (7) 

linearly diverges at p — > oo, which precludes the straightforward evaluation 
of these relativistic corrections. 

The crucial observation is that the true relativistic expression for the 
annihilation kernel does not grow up at p — > oo, as distinct from its expansion 
in p/m. So, the initial integral (|5]) in fact converges. 

When treating relative corrections to the decay amplitude, it is convenient 
to single out from it the factor ip(r = 0) = (7ra 3 ) _1,/2 , and a trivial overall 
dimensional factor from M(p). So, from now on we investigate, instead of 
(HP, the following expression: 

r dp 8%a 3 „ W ^ N 
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with dimensionless M(p). 

Let us consider first an auxiliary integral 



dp 8ira 3 



(2vr) 3 p 4 a 4 



[M(f)-M{0)}, 



(9) 



which converges both at low and high p. After the angular integration, the 
dimensionless kernel M(p) depends on the ratio y 2 = (p/m) 2 only, and the 
expression @ reduces to 



- r%[M{y 2 )-Mm- 
7r jo y l 



(10) 



In other words, this auxiliary integral is of first order in a/ir, and there- 
fore of no interest for our problem. This is a first-order radiative correction 
absorbed already by — 10.28a /tc in (|2|) (for orthopositronium). In fact, we 
have neglected in this argument the kernel dependence on the positronium 
binding energy (it certainly exists at least in the noncovariant perturbation 
theory we were starting from in ||10|| ). But corrections effectively neglected 
in this way, are of higher odd powers in a. 

So, expression (§) can be used gulator, and in the now rapidly 

converging integral 



dp 



(p 2 a 2 



p 4 a 4 



[M (p) - M (0)] 



(11) 



we can safely expand M(p) up to (p/m) 2 included. In this way we obtain 



(2tt); 



(p 2 a 2 + l) 2 p A a A 



~4 a 



;i2) 



This is an alternative derivation of the master formula (|j) used in our article 
[ [Tcfl . In our opinion, this derivation leaves no doubts in the correctness of 
this prescription. 



One of us (I.Kh.) is grateful to P. Labelle for discussions. 
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